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Recently, Seebeck coefficients of ferromagnetic conductors are found to be spin-dependent. How¬ 
ever straightforward method of accurately determining its spin polarization is still to be developed. 

Here, we have derived a linear dependence of anomalous Nernst coefficient on anomalous Hall angle 
with scaling factor related to spin polarization of Seebeck coefficient, which has been experimentally 
verified in [Co/Pt] n superlattices. Based on the dependence, we have also evaluated spin polariza¬ 
tion of Seebeck coefficient of some ferromagnetic conductors. Besides, we have also found a new 
mechanism to generate pure spin current from temperature gradient in ferromagnetic/nonmagnetic 
hybrid system, which could improve efficiency from thermal energy to spin current. 
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Since the birth of spintronics, methods as nonlocal spin 
injection [1] , spin Hall effect [2, 3], spin pumping via fer¬ 
romagnetic resonance [4, 5] and circular polarized optical 
excitation [6] have been developed to generate pure spin 
current (J s ), a long-aspired entity for its versatile func¬ 
tionality to realize magnetization rotation [7, 8], rever¬ 
sal [9-11] or domain wall motion [12]. Recently another 
novel category of methods, spin Seebeck effect (SSE) [13- 
20] and spin-dependent Seebeck effect (SDSE) [14, 21], 
utilizing VT to produce pure J s were demonstrated. VT 
leads to in a ferromagnetic (FM) layer out-of-equilibrium 
distribution of magnons which further damps back to¬ 
ward equilibrium state by emitting J s into an adjacent 
nonmagnetic (NM) layer in SSE [14]. In SDSE [21], in¬ 
stead, VT and resulting Seebeck voltage induced by the 
VT leads to, respectively, diffusion current and drift cur¬ 
rent. They not only have opposite directions but also 
different spin polarizations, thus generating a pure J s 
depending on — S±. Here SMn is Seebeck coefficient 
of spin up (spin down) carriers. In order to evaluate 
S t - S i as well as P s = (S t - Si)/(S t + S 4 ), Slachter [21] 
introduced the J s into a NM material, measured spin 
accumulation in the NM material from which they fur¬ 
ther backward calculated J s . However, transparency 
of J s from FM to NM layers is hard to evaluate accu¬ 
rately due to difficulty in determining spin mixing con¬ 
ductance [22] and some other interference effects such as 
interface spin memory loss effect [23] and magnetic prox¬ 
imity effect [19]. Here, instead, we have adopted another 
strategy, making use of inverse spin Hall effect (ISHE) of 
FM layer itself, to measure the J s inside FM layer. We 
will also show in the following (1) that anomalous Nernst 
effect (ANE) whose scenario has long been in mystery 
has close relation with SDSE and anomalous Hall effect 
(AHE) and (2) that Ps is obtainable by measuring ANE 
and Seebeck coefficient as well as anomalous Hall an¬ 
gle. These findings would be beneficial for developing 
a straightforward method to evaluate Ps and thus also 
for searching materials with large transforming efficiency 


from heat flow to pure spin flow. 
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FIG. 1. Schematic diagrams of (a) measurement geometry 
and (b) physical scenarios of SDSE and ANE. Diffusion cur¬ 
rent is totally balanced by drift current while their spin coun¬ 
terparts is not, leading to a pure spin current along x and 
further a charge current along y due to ISHE in FM. 

[Co(12/n nm)/Pt(12 /n nm)] n with n = 1 — 8 were 
magnetron-sputtered on Si/Si0 2 (500 nm) substrates 
with a base pressure of 1.0 (10 -6 Pa) at room tem¬ 
perature. The number of Co/Pt interfaces increased 
as (2n — 1) while their total thickness was kept as 12 
nm. Co(12 nm), Pt(12 nm) and X(20 nm)/Pt(2 nm) 
(X=NiFe, CoFeB, Co, Fe) were fabricated for compar- 
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ison. Then all films were patterned as Nernst Bar as 
shown in Fig. 1 (a). In Nernst measurement, heating 
current was transported between Pads A and B, result¬ 
ing in a thermal gradient in transverse direction. In this 
direction, no current flowed between Pads C and E since 
the subtrees were isolated to the heating bar by 7 pm. 
Then Pads C and D were used to detect a pure Nernst 
voltage free of voltages due to any current related effect 
such as Peltier Effect or Anomalous Hall Effect. Pads C 
and E were used to detect Seebeck voltages. Magnetic 
field was provided by Physical Properties Measurement 
System (PPMS-9T, Quantum Design). 

Spins in Co were first aligned along z by a large 
B z . W X T was generated by applying I between Pads 
A and B. Seebeck voltage E c was then built. W X V C = 
+ <J^S ^)/(JcoVT = S'coV^T. Here , a± and 
(JCo = + <Ji were conductivity of spin up, spin down 

channels and total conductivity, respectively. Although 
diffusion current crcoSco^xT was fully balanced by drift 
current crcoV^Ec, diffusion spin current Pt^CoPcoV^T 
differed from drift spin current Perco^xVc due to their 
different spin polarizations, thus J s emerging with J s ^ x = 
(P — Pt)&Co^xV c (SI A) (Fig. 1 (b)). Spin polarization 
P of conductivity, spin polarization Pt of thermoelectric 
conductivity (crS) and spin polarization Pg of Seebeck co¬ 
efficient were defined as P = (oq —cyQ/crco, Pt = (cr^S^ — 
cr;P 4 ,)/(cr t P t + and P s = - 5j,)/(5 t + S±), re¬ 

spectively. They correlated via Pg = (Pt—P)/( 1 — PpP). 
The derived J s was the same as derived by Slachter [21]. 
They further introduced J s into a NM material and mea¬ 
sured it via spin accumulation in the NM material while 
we calibrated it by ISHE in FM itself directly (SI A). 
Resulted voltage En (ANE voltage) was thus 

V y E N = (1 - Pt/P)0^ x V c (1) 

ANE coefficient rj = VyVp/(poM 0z V x T), po and M 0 


permeability of vacuum and saturated magnetization, re¬ 
spectively. Anomalous Hall angle Op, ANE coefficient p 
and Seebeck coefficient Sco were correlated via Eq. (2). 


r] = (1 - Pt/P)OuS C o/(^M 0z ) ( 2 ) 

Here Op = PO s where 0 S , spin Hall angle, characterized 
spin-orbit coupling (SOC) strength of FM material. Lin¬ 
ear scalability of p with Op and S was obtained previously 
in nonmagnetic materials [24] but not in ferromagnetic 
counterparts. It deserved special attention that quantity 
(V^Vn/V^T— 0pSc o ) which was lack of deep comprehen¬ 
sion in previous literatures [25, 26] became interpretable 
after Pt was introduced. Similar with Op, Nernst angle 
#n = Pr£ s . 

Multilayer differed subtly from a single layer in ther¬ 
motransport properties. Taking Pt and Co bilayer as 
an example and supposing they had the same thickness, 
Peff = (tfCoSco # crptPpt)/ (p Co + cr Pt ). J S:X in Co of the 
bilayer was remarkably modified as J s , x = aco/ (^Co + 
*Pt)Sco(P-PT)*CoV x T+(PS Pt -P T Sco)(TCoCTPt/((TCo + 
c r Pt)V cc T (SI B). There existed two origins contributing 
to J s?a ,, the 1 st term from intrinsic SDSE but reduced 
by a factor aco/ (ct.Co + crpt) due to shielding effect of Pt 
and the 2 nd term named as extrinsic SDSE from varia¬ 
tion of Seebeck coefficients in adjacent layers. To high¬ 
light the extrinsic one, one could suppose a fictional sys¬ 
tem where intrinsic SDSE vanished (Pt = P). Then 
Js,x = -P(S Co - Spt)<7Co<7pt/(<7Co + crpt)V X T . In con¬ 
trast, J s , x = -2(Sf - P|)oyoq/((7 t + cri)\7 x T for a single 
FM layer. If Co and Pt in the bilayer were deemed as 
two spin channels with P and 0 spin polarizations, re¬ 
spectively, the bilayer behaved similarly to a FM single 
layer. SOC considered, spin current also introduced a 
transverse ANE voltage via Eq. (3). 


V^En _ (1 - Pt/P)^Co/(^Co + cpt) + (Ppt/Pco Pr/P)^Pt/(^Co + apt) 
{0pS7 X V C ) 1 + (crptSpt)/ {cTCoSco) 


The above analysis showed that ANE, AHE and SDSE 
were correlated with each other: SDSE in ferromagnetic 
materials generated pure spin current which led to a 
transverse ANE voltage through SOC in the same man¬ 
ner as occurred in AHE. It also indicated that Pr/P was 
obtainable from anomalous Nernst effect. 

We have prepared pure Co films (12 nm) and [Co(12/n 
nm)/Pt(12 /n nm)] n superlattices to modulate SOC 
strength, then measured their AHE, Seebeck effect (SE) 
and ANE in order to measure Pr/P of pure Co as well 
as to testify existence of the extrinsic SDSE. Practically, 
Seebeck voltage and Nernst voltage were measured be¬ 


tween Pads C-E and Pads C-D, respectively (Fig. 1 (a)). 
To remove a same parameter (ViT)/[(Ti—Xe)/he] intro¬ 
duced by nonlinear V X T (SI C), Eq. (3) would be divided 
by Eq. (1). Here (ViT) is temperature gradient at Point 
I while (Ti — Pe)/he is average temperature gradient be¬ 
tween Point I and Pad E. 

Fig. 2 (a) showed B dependence of p xy of the samples 
(SPn, n = 1 — 8) at 300 K. 2.5 T Field was high enough to 
saturate M along z. p xy and p xx both linearly depended 
on the number of Co/Pt interfaces (2 n — 1) (Fig. 2 (b)), 
indicating interfacial scattering and extrinsic skew scat¬ 
tering gradually dominated p xx and p xy , respectively. Op 
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FIG. 2. (a) B dependence of p xy , (b) p xx and p xy of SPn, (c) 
Oh of SPn at 300 K and (d) scaling law between p xy and p xx 
of SP1 and SP8 from 300 K to 10 K. The black and red lines 
show linear fitting of p xy vs. p xx , respectively. 


of SPn was summarized in Fig. 2 (c). Oh of SP8 was 
0.040, much larger than 0.0075 of SP1 and 0.011 of pure 
Co. Scaling law p xy oc p xx @ was measured by varing T. 
Scaling exponent (3 of SP1 was 2 while that of SP8 was 
reduced to 1.4 (Fig. 2 (d)), also demonstrating dominat¬ 
ing role of skew scattering in samples with larger n [27]. 
Similar results that Co/Pt interfaces induced skew scat¬ 
tering were also reported by Canedy [28]. 

Seebeck coefficients of pure Co (12 nm), Pt (12 nm) 
and the superlattices were also measured with heating 
I flowing through Pads A and B and voltmeter con¬ 
necting Pads C and E (SI D). V X T was proportional 
to heating power P . Thus v c = dV c / d P and z/n = 
dV^/ d P were used to indirectly characterize S and p. 
n c ,p t = -(7.8 ±0.4) qV/W and z/ c?Co = -(23 ±2) qV/W. 
S Pt = -5.0 qV/K [29]. Then S C J= -(14.7±1.5) qV/K, 
close to —22 jllV/K applied in Ref. 29. The z/ c?n exhibited 
no systematic dependence on n. Averagely z/ c?n = (16 ± 
3.0) qV/W. Accordingly, S n = —(10.3±2.0) qV/K. The 
formula S n = {acoSco ± 0‘PtSpt)/(o'Co ± opt) predicted 
S n = —(10.1 ±0.8) qV/K considering (Jpt/(Jco — 1-1, 
consistent with the measured value. T\ — Tp and average 
V X T between Point I and Pad E were thus 0.78 K and 22 
K/cm at P = 0.5 W. Even higher \7 X T existed in regions 
closer to the heating bar (SI C) where giant ANE was ex¬ 
pected. Due to effectiveness to produce large VT, this 
geometry has also been applied to study thermoelectric 
properties of carbon nanotubes [31] and graphene [32]. 

Vn evolved with B z in a similar manner with Eh and 
sign of Vn was reversed as \7 X T was reversed, a typical 
feature of ANE (Fig. 3 (a)). Fig. 3 (a-b) also showed or¬ 
dinary Nernst effect was small enough to be ignored up to 
3 T. Saturated En increased linearly with P as expected 
(Fig. 3 (b)), thus pn = cIEn/cLP did not depended on 



FIG. 3. (a) En between Pads CD and FG, inset showing 

corresponding setups, (b) the dependences of En between CD 
on B under different P, (c) the dependence of z/n and Oh on 
2n — 1 and (d) the scaling relation between nn and Oh , the 
lines representing linear fittings of z'n vs. Oh as n < 5 and 
n > 5. 


P . SP8 exhibited much higher z/n than SP1 at the same 
condition. Similar with linear dependence of Op on n 
as n < 6, z/n of the superlattices also linearly increased 
with n (Fig. 3 (c)). z/n ,5 = 2.6z/n,i, indicating interfacial 
SOC was the common origin of AHE and ANE (Fig. 3 
(d)). This linear relation also validated Eq. (1) and Eq. 
(3). Though dependencies of #h and z/^ on n deviated 
linear relation beyond n > 5, z/^ still linearly scaled 
with #h as n > 5 with a larger dv^/ d$H- du^/ d#H = 
(0.33 ± 0.02) mV/W and (0.80 ± 0.06) mV/W, respec¬ 
tively, as n < 5 and n > 5. 

Co = z/n,Co/(#hz±,Co) = 42 ± 4 for the Co single layer 
while C n = (dz^sr/ d#H)/±c,n = 21 ±4 as n < 5, averagely. 
As indicated by Eq. (3), C n should have been about 
40%C 0 = 17 ± 2 as cr Pt ^ 0.91 <t Co and 5 Pt /5co = 0.34 if 
only the intrinsic SDSE contributed to spin current. The 
measured C n (n < 5) was (1.24 ± 0.28) times as large as 
the estimated value, indicating the intrinsic SDSE played 
major role. In contrast, C n = 50 ± 10 as n > 5, two 
times larger than the expected value. The extra part, 
as per our knowledge, has to be attributed to the ex¬ 
trinsic SDSE. By calculating [Eq. (3)/ Eq. (1)] with 
C n = 50 ± 10 and Co = 42 ± 4, one could remove the 
parameter (ViT)/[(Ti — Te)/t*ie] (SI C) and then ob¬ 
tain P T /P = 1.56 ± 0.35 and P T = -(0.62 ± 0.14) if 
P = —0.4 [33]. According to the results in Ref. 29, 
Pt was estimated as 51%, also close to the result mea¬ 
sured here. According to Eq. (2), p = —50 nV/(K T) 
considering /i 0 M 0 = 1.7 T, in the same order of p of 
Co 2 FeAl [34]. P s = -(0.30 ±0.20) for pure Co and V X T 
at Nernst cross was about 0.165 K/qm at P = 0.5 W 
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and (ViT)/[(Ti — Te^/vie] was about 75. For a pure 
Co layer, dJ s ^ x /d\7 x T = (11 =b 7.5) A/(mK) there¬ 
fore, J s ^ x = (1.8 ± 1.2) MA/m 2 as P = 0.5 W, far 
enough to drive rotational and translational motions of 
skyrmions [35]. It was noteworthy that the extrinsic 
SDSE arose as n > 5 or spacing between Pt layers be¬ 
came shorter than 2.4 nm, hinting the extrinsic SDSE was 
merely remarkable in interfacial regions by less than 1.2 
nm. The length scale is comparable to typical spin dif¬ 
fusion length in Pt. Compared with the intrinsic SDSE, 
the extrinsic SDSE could be more flexibly modulated by 
suitable choice of materials whose Seebeck coefficients 
differed by large values. 
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FIG. 4. (a) AHE and (b, c and d) ANE between Pads F-G 
of X (Co, Co 40 Fe 40 B 20 , and Ni 80 Fe 20 , 20 nm)/Pt(2 nm) at 
elevated P. 


We also investigated Pt /P of some other ferromagnetic 
materials as Co 40 Fe 4 0 B 20 , Ni 80 Fe 20 (Fig. 4) and Fe (not 
shown here). r]F e was negligibly small due to very small 
and positive Ppe- Structures of X (20 nm)/Pt (2 nm) 
were deposited where X=Co, Co 4 oFe 40 B 20 and Ni 80 Fe 20 . 
Pt layer functioned as anti-oxidation layer. The u c 
of Co, Co 40 Fe 40 B 20 or Ni 80 Fe 20 were —19.4 pV/W, 
—7.78 pV/W and —16.1 pV/W, respectively. Their z/ N 
were 7.27 pV/W, 21.7 pV/W and 11.3 pV/W, respec¬ 
tively. Their Hall angles were 0.011, 0.048 and 0.0032, 
respectively. Thus their |Cx/pt| were about 34, 58 and 
220, respectively. (ViT)/[(Ti — Te)/ge] = 75 taken into 
account, their Pr/P were accordingly about 1.45, 1.77 
and 3.93. If |P| of Co, Co 40 Fe 40 B 20 and Ni 80 Fe 20 were 
taken as 0.4 [33], 0.5 [36] and 0.2 [37], their |P T | were 
about 0.58, 0.89 and 0.79, respectively and their P$ were 
about 0.23, 0.70 and 0.70. |Pr,Co/Pt| = 0-58 was compa¬ 
rable to |Pr,Co| = 0.62 db 0.14 measured in the superlat¬ 
tices. Ni 80 Fe 20 and Co 40 Fe 40 B 20 had higher Pp and Ps 
than Co. The result that Ni 80 Fe 20 had larger Ps than 
Co has also been reported in Ref. 29, though different in 
magnitudes. 

In summary, ANE was interpreted as joint action of 


SDSE and SOC. ANE, AHE and SDSE were then com¬ 
bined to characterize Pt and Ps in ferromagnetic films. 
For cobalt, P T = -(0.62 ±0.14) and P s = -(0.30±0.20) 
and dJ S)2C / dV^T = (11±7.5) A/(m K). For Co 40 Fe 40 B 20 
and Ni 80 Fe 20 , Pt/P was about 1.77 and 3.93, respec¬ 
tively. ANE could provide a straightforward method to 
calibrate spin polarizations of Seebeck coefficients of fer¬ 
romagnetic conductors. We also distilled an extrinsic 
SDSE originating from variation of Seebeck coefficients. 
It only significantly affected spin transport within interfa¬ 
cial regions of less than 1.2 nm for Co/Pt interface. This 
extrinsic mechanism could offer another flexible strategy 
to generate pure spin current via VT, comparable to in¬ 
trinsic mechanism, which could be applied especially in 
skyrmion spin-caloritronics due to its moderate efficiency 
from heat to spin current. 
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A. Derivation of relation between ANE co¬ 
efficient, AHE angle and Seebeck coefficient in 
single ferromagnetic layer system 

We will follow the semiotic system in Ref. [SI] in the 
following discussion. According to diffusive and linear 
response theory [SI], charge current density J c , spin cur¬ 
rent J s and heat flow Q are driven by electric potential 
gradient V/x c , spin potential gradient V/x s and tempera¬ 
ture gradient VT, according to Eq. (SI). 


Jc 


' 1 P ST ' 


V/z c /e 

Js 

= a 

P 1 P t ST 


V/U s / 2 e 

Q 


ST P t ST kT/<t 


-VT/T 


Here J c / S = J^- ± where J^ and J ^ are current den¬ 
sity in spin up (t, defined according to direction of mag¬ 
netization) and spin down (j.) channels. S', a and k 
are total Seebeck coefficient, total conductivity and total 
thermal conductivity of a ferromagnetic material, respec¬ 
tively. S = (Sfa+ + S\ (J\ )/(<7t + <J\ ) and a = + cr \. 

P=(a t -a-)/U + a i ). pi = [d (Pa)/ d e U/^/Yd 

oq and oq are spin-dependent conductivities. Electro¬ 
chemical potential fi c = + /xjJ /2 and spin accumu¬ 

lation /x s = (/x-f — yiq )/2 where fi^ and fii are electro¬ 
chemical potential for spin up and spin down electrons. 
Compared with (i c , is relatively small in bulk region as 
implemented too as a hypothesis in Ref. [S2, S3]. Then 
we arrive at Eq. (S 2 ). Q is ignored in Eq. (S 2 ) since heat 
flow driven by VyU c and V// s in our experiment is much 
smaller than what we have applied by an external heat¬ 
ing current of about several tens of milliamps. In other 
words, VT is dominantly determined by external heating 


current instead of by V/x c and V fi s . 


'Jc 


V/ic/e - .SVT 

A 

— < 7 

PV/i c /e - PtS'VT 


(S 2 ) 


For an open-circuited film, J c = 0, thus VyU c = eSVT, 
/x c /e being Seebeck voltage actually in this case and Eq. 
(S3) holds. 


J s = (P - P t )ctSVT = {P- c /e (S3) 


Interestingly, spin current J s still survives even though 
net charge current J c in this case does not exist if 
P 7 ^ Pt, which is the reason why the J s is named as a 
pure spin current. Actually, Pr = (oqSf — oqS 4 ,)/(oqSf + 
oq Si) taking Mott relation = ^-(|)(fcT)^ into ac¬ 
count where k is Boltzmann constant. Pr is spin polar¬ 
ization of thermoelectric conductivity (a = aS). Thus 
Eq. (S3) could also be equivalent to equation (S4), which 
is also obtained in Ref. [S3]. 

Js = -(1 - P 2 )/ 2 (S t - SJaVT (S4) 


The above equation tells us (1) that pure spin cur¬ 
rent originates from difference between S\ and and 
( 2 ) that spin polarization of Seebeck coefficient Ps = 
(S^ — Si)/(S^ + Si) could be obtained if the pure spin 
current J s generated by VT could be measured. The 
relation among Pr, Ps and P is as shown in Eq. (S5). 


Pr 


Ps + P 
1 + PsP 


(S5) 


Pure spin current vanishes only if S^ = Si or Ps = 0 or 
P T = P. 

In this paper, in order to detect J s , we have adopted 
spin-orbit coupling (SOC) in ferromagnetic material 
(FM) itself to transform the J s into a charge current 
J c . Before going deep into anomalous Nernst effect, we 
first discuss how a longitudinal J s induces a transverse 
J c by anomalous Hall effect (AHE). Supposing current 
J is applied in FM along x, the J could be decomposed 
into two spin components: J (1 + P)/2 for spin-up elec¬ 
trons and J(1 — P)/2 for spin-down electrons. Thus a 
spin current PJ is produced. Without loss of any gen¬ 
erality, we further suppose electrons with different spins 
could have (but not necessarily) different spin Hall an¬ 
gles Oyi. Pq = ( 6 ^ — 0i)/{0i + Qi) and let average spin 
Hall angle 0 S = (Pq + Pq)/ 2 . Thus a transverse current 
J c = 0 t J(1 + P)/2 - 0i J (1 - P)/2 = 0 s J{P + P e ) is gen¬ 
erated. Anomalous Hall angle #h = J c /J = 0 S (P + Pq). 

As occurred in AHE, spin current generated by VT 
in spin-dependent Seebeck effect (SDSE) could also be 
transformed as a transverse charge current through SOC, 
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leading to anomalous Nernst effect (ANE). However, dif¬ 
ferent from AHE, J s in SDSE is a pure spin current which 
leads to a transverse current of J c = 0^ J s /2 + 0^ J s /2 = 
0 S J S . Ratio of J c /J s = 0 S in this case, deviating from On 
by a factor of (P + Pq). 

Considering (S3), the induced Nernst voltage (Hn), in 
open circuit condition, is then 

W N = O s J s /a = (P - P T )0 s \7fi c /e (S 6 ) 

It is worth noting that anomalous Hall angle On equals 
to 6 > S (P + Pq). Thus 


VF n = ^feV Mc /e 

-r + re 

c _ VVk = P-Pt 
~ c /e) P + Pe 

Anomalous Nernst voltage Vn proportions product of On 
and V/i c /e with scaling factor C depending on spin polar¬ 
ization of conductivity, thermoelectric conductivity and 
spin Hall angles in FM. Here 0^ = or Pq ~ 0 is sup¬ 
posed because ( 1 ) we think electrons with different spins 
experience the same scattering potential in FM and (2) 
there have been few experimental or theoretical litera¬ 
tures to show their difference according to our knowledge. 
Thus 


W N = E_p£ H v Mc /e 

(S7a) 

W N P ~ Pt 

~ OnVUJe) P 

(S7b) 


The above results could also be derived from another 
formalism as adopted in Ref. [S2, S4]. Supposing = 
cTijEj — OLi k VkT , a and a being conductivity and thermo¬ 
electric conductivity tensors, respectively. Ji and Ej are 
charge current density along i and electric field along j, 
respectively. Under insulative boundary conditions and 
V y T = 0, 


Px 

E- 

$ 

> 

ii 

Pxx^xx T Pxy^xy 

_Ey_ 


Pxy&xx Pxx&xy 



p is resistivity tensor. 

Seebeck tensor is defined as usual as Sij = P^/VjT . 
Then 


Sxx 


Pxx^xx T Pxy^xy 

Syx 


Pxy&xx Pxx&xy 


(S9) 


Hall angle tan#H = Pxy/Pxx • For convenience, we define 
another angle, Nernst angle, as tan#N = OL X y/oL X x • Then 


SxX 

Syx 


1 + tan On tan On 
tan On — tan On 


Bear in mind that On and On are both much smaller 
than 1 for majority ferromagnetic materials, thus 
S xx = Pxx^xx and Sy X = S xx (tan On tan$N). The 
Nernst angle On reflects SOC induced deflection of 
diffusion current under VT. Similar with On = P# s , 
On = P T 0 S . Here, Pp is spin polarization of thermoelec¬ 
tric conductivity defined above. Thus On/O n = Pp/P. 
Then S yx = S xx (0n — On) = S xx 0n( 1 — Pt/P) as shown 
in Eq. (S7). Here we also suppose Pq = 0. 

B. Derivation of relation between ANE coef¬ 
ficient, AHE angle and Seebeck coefficient in 
ferromagnetic/nonmagnetic bilayer system 

Thermal transport properties of ferromagnetic and 
nonmagnetic bilayer system show some subtle differences 
from the case of the single ferromagnetic layer. The bi¬ 
layers are restricted to have the same thickness h and 
the same width w for the ferromagnetic and nonmag¬ 
netic layer for simplicity. In this system, charge and spin 
currents satisfy Eq. (S10). Here subscript Pt and Co 
are used to label quantities in Pt and Co layer respec¬ 
tively. Subscripts f and f also label quantities in Co 
layer. V c = V/i/e Seebeck voltage. 


Jc = — — cr^SfVT — crp t SptVT + (<tc 0 + c r Pt)VU c = 0 

J s ,Co = a t P t VT - <t^VT - (<r t - arf VU C (S10) 


Then effective Seebeck coefficient is modified as 

+ o-ptSpt crcoSco + crptSpt 

^eff = -:-= -:- 

tfCo + crpt cr Co + CTpt 

Similar with the case of single layer in charge trans¬ 
port, the effective Seebeck coefficient is the summation 
of weighted Seebeck coefficient of each layer and each 


spin channel with the weight factor being their conduc¬ 
tivity ratio of However, thermally activated 

spin transport in the bilayer is remarkably different by 
the single layer as shown in the second term of right hand 
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side of Eq. (Sll). 

(P - Pt)(Tc 0 Sco Tjrr , (PS Pt - -PT5'co)(7'Pt(TCo vvrr 

J S) Co = -;-Vi H-;-Vi 

CCo + CPt CTCo + <7Pt 

(SU) 

Another contribution to pure spin current arises due to 
variation of Seebeck coefficient between Co and Pt. To 
highlight its influence, one could image a fictional system 
where intrinsic spin current induced by spin polarization 
of Seebeck coefficient of Co is zero or Pt = P- In this 
system, 

J s ,Co = P(S Pt ~ Sco ) aPtaC ° VT (S12) 

tfCo + cr p t 

Diffusive current ocoScoVT driven by temperature gra¬ 
dient completely cancels drift current crc 0 VW in a single 
magnetic layer. However, the equilibrium between drift 
current and diffusion current inside cobalt layer is broken 
by another adjacent layer of platinum by shielding part 
of drift current inside the cobalt layer. At this time, 
the drift current inside the cobalt layer is (JcoVV^ = 
OCoPeffVT = CTCo(^CoPco + ^Pt^Pt) VT/((TCo + 0Pt) 
while the diffusion current in it still maintains the value 


of acoPcoVT as in the single cobalt layer. There¬ 
fore a net current a Co {^CoSco + crptPpt)VT/(cr Co + 

0pt) - vCoScoVT = (Sp t - S Co )crcoCrptVT/(aco + crpt) 

arises in the Co layer, which is accurately balanced 
by another net current crp t VW — tfptSptVT = (Sco — 
Ppt)c r CoC r PtVT/(crco + <rpt) in the Pt layer for the whole 
bilayer system. However the net current in the Co 
layer is spin polarized while the counterpart in the Pt 
layer is not. Thus a pure spin current of P(5pt — 
5'co)c r CoC r PtVT/(cr Co + <7pt) is generated inside the bi¬ 
layer system as Eq. (S12) predicts. The first term in Eq. 
(S12) originates from the spin-dependent Seebeck coef¬ 
ficient of Co layer. We name this term as the intrinsic 
spin-dependent Seebeck effect.The second term relates 
with variation of normal Seebeck coefficient in different 
layers, which is named as extrinsic spin-dependent See¬ 
beck effect thereafter. 

SOC taken into account, the spin current would also 
induce an anomalous Nernst voltage. Here we have fur¬ 
ther supposed that h is smaller enough than a feature 
length A characterizing length scale where SOC effect of 
Pt has significant effect on spin current in Co, then 


W N 


Ms,Co 
0 Co + 0 Pt 


(1 - Pt/P)^ 0 Sco + (gpt/gco - Pt/P WcoSco 
(c-Co + Opt ) 2 


OrVT 


(S13) 


The relation #h = PO s has been applied. Then one 


could obtain Eq. (S14) considering relation between See¬ 
beck voltage and temperature gradient. 


OuW c 1 + (crptSpt)/ {^coSco) 


It should be pointed out that the above theoretical 
results (S14) are based on the assumption that h is 
smaller than A. If the condition is not satisfied, the 
extrinsic spin-dependent Seebeck effect is supposed to 
disappear. 

C. Estimation of temperature distribution 

Measurement setup is shown in Fig. 1(a) and Fig. 
3(a) inset. Heating current is applied along a bar from 
Pads A to B. Thermal conductance of silicon, [Co/Pt] n 
multilayers and silicon dioxide are estimated in the or¬ 
der of 10 -1 W/K (lateral thermal transport), 10 -6 W/K 
at most (lateral thermal transport) and 1 W/K (verti¬ 
cal thermal transport), respectively, though their ther¬ 
mal conductivities are nearly in the same order or dif¬ 
fer from others only by 1 ~ 2 orders of magnitude, 


150 W/mK for silicon [S5], 70 — 100 W/mK for met¬ 
als [S6, S7] and 1.5 W/mK for Si0 2 [S8, S9]. Large dif¬ 
ference in thermal conductance of different layers results 
from large difference in thickness of each layer. There¬ 
fore heat flow generated by heating bar tends to first pass 
through Si0 2 layer underneath the heating bar and is 
then transported through silicon substrate into environ¬ 
ment. It is mainly silicon substrate that dominates heat 
flow transport and that determines temperature distri¬ 
bution along the [Co/Pt] n superlattices. Another fac¬ 
tor affecting the temperature distribution is geometrical 
dimensions of sample. Core structure for Seebeck and 
Nernst measurement along x is within 200 pm while di¬ 
mensions of silicon wafers are 5 mmx 5 mm x 550 pm, and 
the measurement structure is placed around the middle of 
the wafer. Considering the large difference in dimensions 
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between the measurement structure and the substrate 
wafer, we suppose the measurement structure is placed 
on an infinite large substrate for simplicity in analysis 
of temperature distribution. Ignoring heat radiation and 
at steady state, one can obtain from energy conservation 
law, 

cpnr dr dT = \knrT\r) — kn(r + dr)T'(r + dr)] d t 

1 rri 

cpr— = — k[T'(r ) + rT"(r)\ = 0 

Here c is heat capacity of silicon, T is temperature, k is 
thermal conductivity and p is density. The distance away 
from the heating bar is r, as shown in Fig. SI. Thus at 
steady state, 

T\r) = - (S15) 

r 

where A is a constant and the temperature gradient VT 
is inversely proportional to r, consistent with Ref. [S10]. 
Eq. (S15) demonstrates that VT could be extremely large 
as the point where we measure ANE voltage is close to 
the heating bar (r —>> 0). Being hard to be realized in 
conventional methods, this large VT could induce a large 
ANE signal, which helps to improve signal-to-noise ratios 
in our setup. We set the distance between the point and 
the edge of the heating bar to be 7 pm . 


Thermally 



FIG. SI. Isothermal line in silicon, temperature at r depends 
linearly on 1/r. 

We also conducted finite elementary modeling (FEM) 
for temperature distribution as shown in Fig. S 2 . As 
predicted by Eq. (S15), large VT indeed locates near 
the heating bar with a very fast decay of VT along x. 
Here all the boundaries were set to be thermally isola- 
tive in FEM. The decay of VT along x here is found to 
be even faster than prediction of Eq. (S15). It should 
be addressed that since VT nonlinearly depends on x, 
VT at Point I (VTj) where we measured Nernst voltages 
would be much larger than average temperature gradi¬ 
ent between Point I and Point E ((Tj — Te)/t*ie) between 
which we measured Seebeck voltages. This is the reason 
why we have not directly applied Eq. (7) or Eq. (S14) 


in our setup to obtain the value of Pp. Here, instead, 
we reasonably assume that similar placement of heating 
bar on large silicon substrates under the same heating 
power would result in a similar temperature distribu¬ 
tion near the heating bar. Under this assumption, we 
could thus use [Eq. (S14)/Eq. (S7)] to cancel the same 
geometry-dependent factor (VTi)/[(Ti — Te^/vie] • Be¬ 
sides, we also applied z/n = dV ^/dP and u c = dV c / dP 
instead of dV^/ dVT and dV c / dVT to indirectly charac¬ 
terize Nernst and Seebeck coefficient since the concrete 
value of VT is hard to be obtained while P could be 
accurately controlled. 

It is worth reaccentuating that large VT, testified by 
large Nernst voltage about 20 pV in our measurement 
and hard to be realized in conventional heating setup, 
could be even further larger in nanodevices where 
spacing between heating bar and detecting circuit is in 
the range of several tens of nanometers in nowadays 
electronics and thus could contribute a higher signal 
which has easier detectability. 



FIG. S2. Temperature distribution of silicon on which there 
exists a heating bar 

D. Dependence of Seebeck voltages on heat¬ 
ing power 

Seebeck voltages of different samples linearly depended 
on heating power (P) as expected since VT linearly 
depends on P. Fig. S3 offset Seebeck voltages of different 
samples for clarity. Their slopes of u c = dU c / dP are 
shown as inset in Fig. S3. Among these samples, Co has 
largest \v c \ about (23 ± 2 ) pV/W while Pt has smallest 
one about (7.8 ± 0.4) pV/W. The \is c \ of the superlat¬ 
tices shows no systematic dependence on the number of 
interface. Their average \ v c \ was about (16T3.0) pV/W. 


E. Magnetic and transport properties of [Co/Pt] n 
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the samples are smaller than 1%, which could introduce 
about 1% at most inaccuracy in heating power by I 2 R. 
Here I is heating current and R is resistance of the main 
bar in Fig. 1(a). It should be pointed out that the heat¬ 
ing power at +3 T is the same as that at — 3 T because 
R is an even function of B at large field. Therefore inac¬ 
curacy in final Nernst voltage introduced by the AMR of 
the heating bar should be even much smaller than 1% and 
thus this AMR effect was ignored in our analysis in the 
main text. The evolution rule of the AMR with increase 
in period number n is beyond scope of this manuscript 
and will be discussed elsewhere. The magnetizations of 
the superlattices are weakly dependent on the period n 
and their average magnetizations are about 1370 emu/cc. 


FIG. S3. Dependence of offset Seebeck voltage V c of super¬ 
lattices on heating power 



Bz (T) 


FIG. S4. AMR of pure Co and superlattices under magnetic 
field along z direction 


multilayers and pure Co layer 

Anisotropic magnetoresistance (AMR) of the superlat¬ 
tices and the Co layer are shown in Fig. S4. During 
AMR measurement, magnetic field is applied normal to 
plane, the same direction as we conduct Nernst measure¬ 
ment and anomalous Hall measurement. The AMR of all 
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